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I. INTRODUCTION 


The relativistic binary collision of two heavy ions in the Relativistic Heavy Ion Collider 


(RHIC) [l] or Large Hadron Collider (LHC) |2j has been significant for understanding of dy¬ 


namics of quark gluon plasma (QGP) jsi]. In particular, the kinetic study of QGP has been 
done to describe the dynamics of QGP in the framework of the relativistic Boltzmann equa¬ 
tion ( 4 ], whereas the perfect fluidity is explained in framework of the relativistic hydrodynam¬ 
ics by Hirano and Gyulassy s|. I 11 a similar way, inelastic relativistic-collisions are described 
by the inelastic relativistic-Boltzmann equation. Actually, the inelastic relativistic-collision 
between two partons is significant for understanding of QGP 6]. Meanwhile, the character¬ 
istics of the inelastic relativistic-Boltzmann equation have not been discussed. In this paper, 
we consider the characteristics of one dimensional (ID) inelastic relativistic-collisions using 
ID inelastic relativistic-Boltzmann equation. We remind that ID elastic collisions never 
yield changes of states of partons, because the binary elastic collision means the exchange 
of the velocity of colliding partons. We, however, consider that ID inelastic relativistic- 
collisions are significant as the first step to understanding of inelastic relativistic-collisions, 


as ID inelastic collisions have been discussed for the granular gas [ 7 ]. The primary aim of 
this paper is an investigation of relativistic effects on time evolutions of the temperature 
and flow velocity under the spatially homogeneous state, when the collision is described by 
the inelastic relativistic-collision. Here, readers must notice that the spatially homogeneous 
state corresponds to the spatially uniform state of partons. Thus, the spatially homoge¬ 
neous state allows the non-zero flow velocity. The main reason why we consider ID inelastic 
relativistic-collision is that analytical result on the time evolution of the temperature can be 
obtained using ID inelastic relativistic-Boltzmann equation, whereas the calculation of the 
time evolution of the temperature via 3D inelastic relativistic-Boltzmann equation involves 
mathematical difficulties. The analytical result of the time evolution of the temperature is 
compared with the numerical result, which is obtained by solving ID inelastic relativistic- 
Boltzmann equation on the basis of the direct simulation Monte Carlo (DSMC) method 8] 
9j. Throughout this paper, we assume that partons are composed of hard spheres with 
equal mass and diameter. Meanwhile, we remind that the collisional differential cross sec¬ 
tion depends on the momentum transferred between two colliding quarks, and the collisional 
deflection-angle also depends on the momentum, which is transferred between two colliding 
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partons 


B, 


id [3l, in the energy-regime of the asymptotic freedom. Additionally, we consider 


the reaction (inelastic relativistic-collision), qq —>■ qq + g (q: quark, g\ gluon), under the 
thermally relativistic state [9], in which we can assume the asymptotic of freedom, when 
the thermal energy of quarks is larger than A QCD 10], namely, thermally relativistic. We, 


however, neglect effects of the spin of quarks to simplify our discussions. Finally, Grad’s 14 
moments {4] are calculated using Eckart’s decomposition 


0 . 


II. TIME EVOLUTION OF TEMPERATURE UNDER EQUILIBRIUM STATE 


Firstly, we assume that the inelastic relativistic-collision qq —>■ qq + g does not have its 
reverse reaction, namely, qq + g —» qq, whereas we can assume another path to the energy 


loss by the radiation in qq collisions 


12] • Here, the successive inelastic relativistic-collisions, 


qq —$■ qq + g, occur, exclusively. In other words, the reactive collision cross section of 
qq + g —* qq is assumed to be zero. Finally, such successive inelastic relativistic-collisions 
are expressed with the inelastic relativistic-Boltzmann equation. 

ID inelastic relativistic-Boltzmann equation is formulated under the spatially homogeneous 
state as 

„0 d f (t>P) 


P 


dx° 


= A 


J 


f (t, p") f (t, p") - / (t, p) f (t, p*) 


, dp* 
P 2 


( 1 ) 


where / (t,p) is the distribution function, in which t is the time, namely, x° = t, p° = mcj(v) 
(7(h) = l/\/l — v 2 \ Lorentz factor, v: velocity of a quark, v = v/c, c: speed of light, m: 
mass of a quark), and p = p 1 = mv^{y) ( p a : two momentum (a = 0,1)). I11 Eq. (1), J 
is the Jacobian owing to the inelastic relativistic-collision. The right hand side of Eq. (1) 
corresponds to inelastic relativistic-collisions, where F = g 0 / (p°p°) (g 0 : Mpller’s relative 
velocity) 13]. A in Eq. (1) defines the total collision cross section. In later discussion, we 
set m = 1 and c = 1 to simplify our discussion (i.e., v — v). 

As a result of the direct inelastic relativistic-collision, momentums of two colliding quarks, 
namely, p and p*, change to p' and p*, which are defined by 


/ 1 1 + A 1 \ 

P =P + — (P*-P) 

P* = P* ~ ^ ip* ~ P) , 


( 2 ) 

( 3 ) 
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where A is the inelasticity coefficient (0 < A < 1). On the other hand, momentums of two 
colliding quarks, namely, p" and p", change to p and p*. in which p" and p" are defined by 


" i 

p = p + 
P* = P* ~ 


1 + A 
2A 

1 + A 
2A 


(P* 

(P* 


P ), 


P) ■ 


(4) 


Consequently, the total momentum is conserved by the inelastic relativistic-collision, whereas 
the total energy (E + E* = \Jl + p 2 + \/l + pi) is not conserved by the inelastic relativistic- 
collision. Finally, J in Eq. (1) is the Jacobian, which is defined by 


J = |det d(p",p")/d(p,p*) | 1 = A. 


(5) 


We fix A = 0 in later numerical analysis. 

The time evolution of N a (= J^° oo p a fdp/p 0 ) is obtained by multiplying p“/p° by both sides 
of Eq. (1) and integrating over the momentum space dp/p° as 


d t N a 



(P a ' P*' P° 



f (t,p) f (t,p*) F 


dp * dp 
p o p o' 


( 6 ) 


In the right hand side of Eq. (6), (;p a '/p 0 ' + p*'/p°J - p a /p° - p“/pl) f (t, p) f (t, p*) F 

0, when a = 0. Here, we must remind that the number density (n) temporally changes in 
accordance with the change of the flow velocity (u) owing to d t N° = d t (nU°) = 0, that is 
defined by two velocity, U a = 7 (u) (1, w). On the other hand, we obtain dtN 1 ^ 0, because 
we obtain p 1 /p 0 ' + p* /p° — p 1 /p° — pl/p® = u' + u'* — u — u* 7 ^ 0 from Eqs. (2) and (3). 

The time evolution of the energy-momentum tensor, T a/3 (= j^° oo p a p^fdp/ p°), is obtained 
by multiplying p a p l3 /p° by both sides of Eq. (1) and integrating over the momentum space 
dp/p° as 


d t T a ? 



( P a, PP*'P* 
\ P 0 ' P* 


pOpP 


pV 


p*p% \ 

p 2 J 


f (t,p) f (t,p*) F 


dp * dp 
p o p o- 


(7) 


In the right hand side of Eq. (7), 

(A/2) [p a y'/p 0 ' + p?'pf/p*' - p 0 p p /p° - p*p£/p*) f ( t,p ) / (f,p*) F&jffy is equal 

to 0 from Eqs. (2) and (3), when a = 0 and [5 — 1. 

Consequently, the time evolution of T 01 is obtained using Eq. (29), when / = J'mj = 
n/(2Ki(x)) exp(— XP a U a ) (x = k6/mc 2 : thermally relativistic measure, k: Boltzmann con¬ 
stant, 6: temperature, K n : ri-th order modified Bessel function of the second kind, Jmj'- one 












5 


dimensional Maxwell-Jiittner function 


Q), 


as 


d t T'l‘ = d, | n¥Mu°U' 


IMx) 


= nU°d, I U 


„ Mx) 

KM 

where we used the relation d t N° = dt (nil 0 ) = 0. 

From Eq. (8), we obtain 

U\t)= c KlMt)) 


= 0 , 


( 8 ) 


(9) 


Mx(t)Y 

where C = U\ 0) £ffjj - 

K\(x(t))/K 2 (x(t)) decreases, as x decreases. As a result, |it| increases, as x ($) increases 
(decreases), whereas |w| decreases, as x ($) decreases (increases). Similarly, x ($) increases 
(decreases), as |w| increases, whereas x ($) decreases (increases), as |u| decreases. Of course, 
u(t) = 0, when m(0) = 0. 


Remark 2.1 


The flow velocity u never be conserved by inelastic relativistic-collisions, even when the 
total momentum of binary colliding quarks is conserved. Provided that f = f mj, M is 
inversely proportion to the temperature and its vice versa, liup^^ U 1 (t) = C indicates that 
U 1 converges to not U 1 ( 0) but C, because lim^oo K\(x(t)) / K 2 (x(t)) = 1 is obtained by 
lim^oo x(t) = oo owing to inelastic relativistic-collisions. 


Next, we consider the cooling process by inelastic relativistic-collisions. Here, we con¬ 
sider the time evolution of T° Q/3 instead of that of T Q/3 to remove terms divided by p° in Eq. 
(7). Additionally, we assume that the distribution function is expressed by Maxwell-Jiittner 
function, namely, / = /mj- 

The time evolution of T^ ij = ffl oo p°p a p^ f m jdp / p° is written as 


d t T 0a P = ^ 


A 


r> OO pOO 


’ —oo J —oo 


p a, p >6 ' + p“Vf - p a pP - p“pP) f M j ( t,p ) /mj (t,p*) 


p o p o 


We introduce two vectors P° and Q a , which are dehned by 


P a = p a + pi 


p a, = p a, +p:', 


Q a = p a -p :, Q a ' = p a> -p*', 


(ii) 
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where we remind that p 1 + p\ = p 1 ' + pi' and p° + p* > p 0 ' + p®' . 

We obtain following relations from Eq. (11) 

P a Q a = 0, P a 'Qj = 0, (12) 

P 2 = P“P q = 4 + Q"Q q =4 + Q 2 . (13) 

From inverse transformation of Eq. (11), we obtain 

p a = \ (P° + Q°), p a . = \ (P“ - Q°), 
i = l (p"‘ + «”'), p:' = \ (p* - Q“') , (14) 

Substituting Eq. (14) into Eq. (10), we obtain 


d t T° E aP = 



,atpp _ pap0^ + 


Q a Q p y\ f mj (t,p) f mj (■ t,p *) g 0 dPdQ,(15) 


where we used the relation Fdp/p 0 dp*/p 0 * = l/2g 0 d,PdQ [13j |. 

The integration of A/S (p a 'P? - P a P^ f M j ( t,p) f M j (t,p *) g 0 dPdQ in the right 

hand side of Eq. (15) is markedly difficult, whereas P 1 ' P 1 — P 1 P 1 = 0. 

Hereafter, we can neglect the integration A/S (p a 'pF — p°pP^j f MJ (t,p) f M j ( t,p *) g 0 dPdQ 

with a — 0 or /3 — 0 in the right hand side of Eq. (15), because we will investigate the time 
evolution of T 011 in later discussion. 

The center of mass system yields relations using Eqs. (2) and (3) 


P a =(P°, o), Q“ = (0,Q), Q al = —A (0, Q ). (16) 


From Eq. (16), we obtain 


Q a 'Q p ' - Q a Q 0 = -Q 2 (l - A 2 ) 



where g a/3 = diag(l, —1, —1, —1). 

Mqllcr’s relative velocity in the center of mass system is 



9 & 


2 


Q_ 

pc 


( 17 ) 


(18) 
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From Eqs. (17) and (18), we rewrite Eq. (15), when a = f3 = 1, as 


d t T° E u 



(V'Q 1 ' - Q 1 Q 4 ) /mj (t, p ) /mj (t p*) g 0 dPdQ 

a r°° r°° / pip i \ jp 

p2-^ J Q'fMJ (t,P) fMJ ( t,p *) 



16Xr(x) 2 V4 + Q 2 
A 


- v n Z* ) Q d dQ 


4Ad(x) 2 7 2 

A 



k 2 (xQ*) UWQ* - ( x g*) Q* ) (Q* 2 - 4) dg* 


n 



K 2 (x) UX 1 - - ^!AA - „n K ( Q.) £) (X! _ 4 ) dx 
X X x) Vx 


= — (l — A 2 ) 

1 ' 4 X tf 1 ( X ) 2 j 2x 

= - (1 - A 2 ) n 2 ^ * [2 X K 2 (2 X ) U'u 1 - g 11 (2 X K 0 (2 X ) + I<i(2 X ))] , 

where Z* and Z* a ^ are defined in Eqs. (33) and (34), Q* = \J Q 2 + 4 and x = Q*/ X - 
From Eqs. (8) and (19), we obtain the time evolution of T E n as 


(19) 


dtlf 1 = - (1 - A 2 ) n 
= — (l — A 2 ) n 


A 


X 3 Ki( X ) 2 

A 

X 3 Afi(x) 2 


2 X K 2 (2 X ) (U 1 ) 2 - r] 11 (2 X K 0 (2 X ) + Ki(2 X )) 


2 X K 2 (2 X )[c 


Kiix) 


V Mx) 


rf 11 (2 X K 0 (2 X ) + Ktfx)) 


20 ) 


We obtain following relation using Eqs. (32) and (8) 

r j pon __ r ± _ 

‘ E v/l + M,MW 

[ 2 X 2 (c 2 x‘ + 8 CV + 3 x 2 + 8 ) KJxYMx? + 2 x ( 6 x 2 - C \ 2 (x 2 - 8 ) + 24 ) A',(. v ) 4 Mx) 

-2 (CY (x 2 + 6 ) - 4 x 2 - 16 ) A,(x) s + X 3 ( 4 C 2 x 2 + X 2 - 8 ) /fi(x) 2 *'o(x ) 3 

-x 5 ^o(x) 5 - ex-AUxJA'otx) 4 ] (x 4 Ai(x) 2 a 2 (x) 3 ) -1 Ad, ( 21 ) 


where we used d t (nil 0 ) = 0. 

From Eqs. (20) and (21), we obtain 


dtX 1 = ~nA (1 - A 2 ) Vh (x, C) X \ 

i’l (x,C) = (1 + C 2 Ab(x) 2 /A^(x) 2 )” 1 [ X 2 tf 2 ( X ) (2C 2 xA' 2 (2x)Ad(x) 2 + (2 X K 0 (2 X ) + K,(2 X ))K 2 ( X ) 2 )] 
[ 2 X 2 (C 2 x 4 + 8C 2 x 2 + 3x 2 + 8) Afi(x) 3 W 0 (x) 2 + 2x (6x 2 - X 4 C 2 + 8 X 2 C 2 + 24) AT 1 ( X ) 4 J A 0 (x) 
-2 (C 2 x 6 + 6C 2 x 4 - 4 X 2 - 16) Ab(x) 5 + X 3 (4C 2 X 2 + X 2 - 8) 

-X 5 Ao(x) 5 - 6 x 4 A' 1 (x)A' 0 (x) 4 ]" 1 • (22) 





















From Eq. (22), we find that the cooling rate parameter 0! (y, C) depends on the frame via 

«( 0 ). 

The left frame of Fig. 1 shows 0i(x, C) versus y and C. 0i(x, C) is symmetric at both 
sides of C = 0. 0i (y, C) approximates to 0, when y approximates to oo, whereas 0 1 (y, C) 
approximates to 1/8, when x approximates to 0. Additionally, 0i (y,C) has a peak, which 
moves toward to y = 0, as \C\ increases. 

Provided that C = 0 (w(0) = 0), we obtain 0! (y, 0) from Eq. (22) as 

2yA0 (2y) + K, (2y) 


01 (x, 0) = 


(23) 


X ( X 2 + 4 ) (x) 2 - X 3 Ko (x) 2 ' 

The right frame of Fig. 1 shows 0i (x, 0) versus y. 0i (y, 0) has its peak value 0.295 
at y ~ 4.14, whereas we obtain linx^o 0i(y, 0) = 1/8 and lim x _ j>00 0x(y, 0) = 2/y / 7fy, 
which coincides with the cooling rate parameter, which is calculated using ID inelastic 
nonrelativistic-Boltzmann equation. 



0 


FIG. 1: 0i (y,C) versus y and C (left frame), and 0i (y, 0) versus y (right frame). 


III. NUMERICAL ANALYSIS OF COOLING PROCESS 

We investigate the characteristics of the cooling process, which is derived from the in¬ 
elastic relativistic-collision between two quarks, numerically. 

Firstly, we set A = 0 and A = 1 in Eq. (1). From Eq. (22), we know that the cooling rate 
depends on y and C, when / = f MJ . As initial data, we consider two tests. One is Test 1, 
in which / is uniformly populated in the range of 0.99 < |u| < 1 at t = 0. The other is Test 
2, in which / is uniformly populated in the range of 0.99 < v < 1 and — 1 < v < —0.9 at 
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t — 0 . As a result, initial distribution functions in Tests 1 and 2 are markedly nonequilib¬ 
rium. Meanwhile, DSMC sampling of the initially equilibrium state / = /mj is markedly 
difficult, when u ^ 0 and ID elastic relativistic-collisions occur. Thus, we use the initially 
nonequilibrium state. 

Figure 2 shows time evolutions of y, which are obtained in Tests 1 and 2. y in Test 1 is 
smaller than y in Test 2 at t = 0, whereas y in Tests 1 and 2 increase with similar in¬ 
clinations in the range of 0 < t < 6. The inclination of y in Test 2 decreases at t ~ 6. 
Provided that / ~ /mj at 6 < t, /q (y,C) in Test 2 is smaller than /q (y, C) in Test 1. The 
cooling rate parameter /q (y,C) increases, as \C\ decreases, when 1 < y, as shown in Fig. 1. 
Provided that \C\ in Test 1 is smaller than \C\ in Test 2, ipi (y,C) in Test 1 is larger than 
ipi (y, C) in Test 2, when / = /mj and 1 < y. In later discussion, we, however, find that 
C = U 1 K 2 (y) /K i (y) is not the temporal constant owing to f ^ fhi.j, because T 01 (^ T^ 1 ) 
includes effects of nonequilibrium moments such as II, q a , and If a/3 |11 |. 



t 


FIG. 2: time evolutions of y in Tests 1 and 2. 

Figure 3 shows / and /mj versus u at t = 17.7 in Test 1 (left frame) and / and /mj versus 
u at t — 177 in Test 2 (right frame). The left frame of Fig. 3 shows that / has higher tails 
than /mj at 0.575 < |u|. Such higher tails are obtained under the homogeneous cooling 
state of the granular gas Q. Meanwhile, / has a higher tail at 0.99 < v < 1 and lower 
tail at v < 0.965. Therefore, the overpopulation of / versus /mj at the high velocity tail in 
both sides of v = u, which is always obtained by the inelastic nonrelativistic-collisions for 
any value of u, is not obtained, when u is similar to the speed of light, as shown in the right 
frame of Fig. 3. The initial marked differences between / and /mj are reduced by inelastic 
relativistic-collisions in both Tests 1 and 2. 
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Symbol: f(v) 
Solid line:/ M /y) 



Symbol: f(v) 



FIG. 3: f[v) and /mj(b) versus v at t = 17.7 in Test 1 (left frame). f(m) and versus u at t = 177 

in Test 2 (right frame). 


Next, we investigate the cooling process of 9, which is analytically obtained in Eq. (22). 
Meanwhile, the calculation of the general solution of y in Eq. (22) is markedly difficult. 
Then, we calculate the limiting solution of y in Eq. (22). Provided that C = 0 [u t=0 = 0) 
in Eq. (22), we obtain solutions of y under two limiting cases, namely, y —> 0 and y —> oo 
from Eq. (23) as 


X{t) = X(0) exp 




(x^0). 

2 

(x oo). 


(24) 

(25) 


Equations. (24) and (25) are plotted together with y in Test 1 in Fig. 4, where A = 1 is 
used in Eqs. (24) and (25). Xt=o hi Eq. (24) is defined by Xt=o hi Test 1 and Xt=o hi Eq. 
(25) is defined by Xt=i 5 in Test 1. Eq. (24) gives a good agreement with y in Test 1 in 
the range of 0 < t < 0.5, whereas Eq. (25) gives a good agreement with y in Test 1 in the 
range of 15 < t. Such good agreements are obtained, because initial nonequilibrium of / is 
reduced by the inelastic relativistic-collisions, as shown in the left frame of Fig. 6. 

Next, we investigate the time evolution of the flow velocity u in Tests 1 and 2. Fig. 5 shows 
time evolutions of u in Tests 1 and 2. u in Test 1 temporally decreases, whereas u in Test 2 
temporally increases, as shown in the left frame of Fig. 5. 
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t 

FIG. 4: Plots of Eqs. (24) and (25) together with time evolution of x in Test 1 . 

As described in Remark 2.1, |w| increases, as x increases. Here, we must remind that 
Ut =o = —e (0 < e -C 1) in Test 1, which is caused by the thermal fluctuation at t — 0, 
yields the time evolution of u in the range of u t=0 < 0. Therefore, the signature of u is 
determined by the signature of u t=0 . u in Test 2 approximates to 0.97, as t increases. Such 
a increase of u via the cooling process of 6 indicates that |u| of quarks are accelerated by 
inelastic relativistic-collisions between two quarks, unless the initial flow velocity Ut =o is 
zero, namely, u t = o = 0 and / = Jmj for 0 < t. In particular, the initial signature of u t = 0 , 
namely, 0 < u t = o or u t =o < 0 coincides with the final signature of lirn^oo u. Of course, such 
a acceleration of |u| is not obtained by inelastic nonrelativistic-collisions, because u is always 
conserved at 0 < t. 



FIG. 5: Time evolutions of u in Tests 1 and 2. 
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Finally, we must mention to effects of nonequilibrium of /, briefly. Figure 6 shows time 
evolutions of II, q 1 , lb 11 ) and U 1 Ki(x) / K 2 (x )? which must be a temporally constant, namely, 
C, when / = Jmj, in Test 1 (left frame) and Test 2 (right frame). As shown in the left frame 
of Fig. 6, U'KM/Kzix) increases in the range of 0 < t < 12 owing to nonequilibrium 
effects and approximates to the constant value in the range of 12 < t. Additionally, the left 
frame of Fig. 6 indicates that II, q 1 and lb 11 ) are significant in the range of 0 < t < 12 and 
temporally damped to zero in the range of 12 < t with similar damping rates. As shown 
in the right frame of Fig. 6, U 1 Ki(x)/K 2 (x) increases in the range of 0 < t < 2.5 owing 
to nonequilibrium effects, decreases in the range of 2.5 < t < 8.8, and approximates to the 
constant value in the range of 8.8 < t. Additionally, II and q 1 are temporally damped to 0 
in the range of 6 < t with similar damping rates, whereas lb 11 ) is temporally damped with 
a slower damping rate than damping rates of II and q 1 . 

Provided that we can assume / ~ f MJ at 6 <t in Tests 1 and 2, jf/ 1 /t 1 (y)/A' 2 (y)| T e a t i •C 
\U 1 K 1 {x)/K 2 {x)\i*t 2 , which is obtained from the left and right frames of Fig. 6, indicates 
that |C| Test ! < |C| Teat 2 might prove (V’i)^ 2 < (V’i)i^t 1 at 6 < t in Fig. 2 by setting t = 6 
to t — 0. 


Test 1 Test 2 



FIG. 6: time evolutions of II, q 1 , II 11 , and U 1 Ffi(x)/ K 2 (x) in Test 1 (left frame) and Test 2 (right frame). 


IV. CONCLUDING REMARKS 

In this paper, we investigated time evolutions of the temperature and flow velocity, when 
the binary collision between two quarks is described by ID inelastic relativistic-Boltzmann 
equation. The absolute value of the flow velocity increases by the decrease in the temperature 



















13 


owing to the inelastic relativistic-collisions, whereas the signature of the flow velocity at 
0 < t coincides with that of the initial flow velocity. The cooling rate of the temperature 


depends on both x (#) and C in Eq. (9), which is the function of the initial temperature and 
flow velocity, under the equilibrium state. To confirm the analytical results, we solved ID 
inelastic relativistic-Boltzmann equation under the spatially homogeneous state. Numerical 


results confirm that the time evolution of the temperature, which is analytically obtained 
in Eqs. (24) or (25), gives good agreements with numerical results of the time evolution 
of the temperature. Of course, the cooling rate depends on nonequilibrium moments such 
as the dynamics pressure, pressure deviator, heat flux and so on, which are nonzero values 
during calculations, as shown in Fig. 6. Additionally, the numerical results surely confirm 
that the flow velocity increases by the decrease in the temperature owing to the inelastic 
relativistic-collisions. Finally, the distribution function does not show the overpopulation 
in the high velocity tail at both sides of the flow velocity (i.e., \u\ <C |u|), when the flow 
velocity approximates to the speed of light. Such a tendency of the distribution function is 
not obtained under the nonrelativistic limit, because the distribution function always shows 
the overpopulation in the high velocity tail at both sides of the flow velocity (i.e., \u\ <C |u|) 
under the spatially homogeneous state. 


Appendix A: Definitions of equilibrium moments 


In this appendix, some equilibrium moments are defined. 
At first, the zeroth order moment is defined as 



(Al) 


In Lorentz rest frame, we obtain 



(A2) 
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Z a ^ 1 2 3 4 5 6 - = j™ oo p a p l 3 p 1 p s ...exp(—xp a U a )^ is obtained by the successive differentiation of Z 
using — ~X.U a as 


= 2K 1 ( X )U a , 

Z af} = 2K 2 ( X )U a U 0 - 2 

X 

Z a ^ = 2K z (x)U a U p U 1 - 2 V + + rf^U a ) ^M, 

X 

Z a W = 2K A (x)U a U p WU 5 

- — a(y) ^ a p wu s + ^jj^s + jjajjS + ^xpxjP + v ^u p U a + rf 0 U a W) 
X 

+ 2 K i x) (rfV 5 + v ai v^ s + rTV 7 * ), 


(A 3 ) 

(A 4 ) 

(AS) 


(A 6 ) 


Similarly, Z* a01& " = f^° oo P a pPp i P s ...exp(—xP 0 ‘U a )j^ is obtained by the successive dif¬ 
ferentiation of Z * using —xJJol as 


Z* = 2 K 0 (Q* X ) , (A 7 ) 

Z* Q = 2 Q*K 1 (Q* x )U a , (A 8 ) 

Z* ap = 2 Q* 2 K 2 (Q*x)U a UP - 2 Q* v af} Kl ^ X \ (A 9 ) 

X 
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